ATOMIC ANTENNA MECHANISM IN HHG AND ATI 



M. YU. KUCHIEV^"), V. N. OSTROVSKY^'') 



(a) ^]^g University of New South Wales, Sydney, Australia 
Institute of Physics, The University of St Petersburg, Russia 



1 Introduction 

This paper reviews recent development of the atomic antenna, a theoretical framework 
which describes a number of laser-induced multiphoton phenomena in atoms. The 
localization of atomic electrons inside an atom drastically suppresses their interaction 
with a laser field. For many processes this circumstance favors multistep mechanisms 
when at first one of atomic electrons is released from an atom by the field. After that 
an interaction of the ejected electron with the laser field results in absorption of energy 
from the field and its accumulation in the form of the electron wiggling energy and ATI 
energy. It is very essential that the energy absorbed by the electron can be transferred 
to the parent atomic core via an inelastic collision of the primarily ejected electron with 
the atom. The collision may trigger a number of phenomena including high harmonic 
generation (HHG), enhancement of above-threshold ionization (ATI), production of 
multiply charged ions. In this physical picture the absorption of energy from the field 
takes place in the region of large separations from an atom, where the electron-laser 
interaction dominates over the electron-core potential. This circumstance results in 
dramatic enhancement of the probability of multiphoton processes. Such a scenario of 
photoabsorption was suggested, apparently for the first time, in Ref. Later the idea 
was rediscovered by several authors in different contexts 0, 0, ^. In current literature 
the above described sequence of events is often referred to as the rescattering, the three- 
step mechanism, or even the simpleman model. The term atomic antenna, suggested 
in Ref. [Q, refers to the fact that the firstly emitted electron plays a role similar to an 
aerial in conventional radio devices, enhancing the absorption. 

It is very important that the physical picture drawn above can be implemented not 
only as a model, but also as a clear and rigorous quantum formalism. In this paper we 
outline two convenient ways to implement the atomic antenna idea. The one originating 
from Ref. is called the factorization technique. It was formulated in detail in Ref. |^ 



and recently applied to HG in Refs. 0. We discuss also another, complimentary 
technique, which is referred to below as the method of effective ATI channels. It is 
close in spirit to the approach developed previously by Lewenstein et a/ P, Q and can 
in turn be linked to the Corcum model 0. From the first glance, these two schemes 
differ very significantly. However, we prove their identity by demonstrating that they 
describe the same physical idea from different points of view. 

The paper is organized as follows. Section ^ describes the factorization technique 
which allows one to present the amplitude of the "complicated" multiphoton process as 
a product of the amplitudes of much more simple, "elementary" processes. Section |^ 
is devoted to the concept of effective ATI channels, which provides important insights 
into the physical nature of complicated multiphoton processes. We derive the effective 
ATI channels from the factorization technique revealing close links between the two 
approaches. Section ^ describes two important examples of "elementary" processes: 
one of them is the photoionization, another one is the electron-atom collision in a laser 
field which results in the generation of the high-energy quanta. These two "elementary" 
processes are vital for description of HG in the framework of the factorization technique. 
Sections |^, || are devoted to several examples illustrating numerically an accuracy of 
methods developed for HHG and ATI. A number of results reported in this paper are 
derived neglecting the Coulomb field of the residual atomic particle which influence the 
active atomic electron. This approach is well justified for negative ions, but needs to 
be modified for the processes with neutral atoms. Section ^ exposes our recent progress 
based on the eikonal approach which allows us to take into account the Coulomb field. 
The concluding Section summarizes the results. The atomic units are used throughout 
the paper unless indicated otherwise. 

2 Factorization technique 

The factorization technique of Ref. |^ can be applied to a number of multiphoton 
processes. In this section our attention is restricted to an important example of HG 
which has been recently considered by the present authors 0. We concentrate on 
the linearly polarized laser field 

F(t) = Fcoscut (1) 

that creates the external potential V{t) = F{t) ■ r acting on atomic electrons. Using 
the second order time- dependent perturbation theory one can present the amplitude of 
HHG (i^ in the single-active-electron approximation in the following form 

d% = -^ Tdt f dt' iUt) I exp(int) U G{t, t') V{t') I Ut')) . (2) 
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Here the brackets ( | | ) imply integration over coordinates, T = 2tt/ uj is the laser period. 
The initial-state wave function is 

0a(r,t) = 0,(r)exp(-zEat) (3) 

where Ea = — is the bound state energy and 0a(r) is the corresponding atomic 
eigenf unction. The potential U = e-r accounts for production of the high harmonic with 



the frequency Q = Nuo and the polarization e. The Green function G{t^ t') describes the 
electron propagation in the intermediate state. The amplitude (Q) takes into account 
only a sequence of events in which HHG follows the absorption of the necessary energy 
from the laser field. A number of omitted, the so called time-reversed sequences, is 
strongly suppressed due to multiphoton nature of the process. 

Neglecting the potential of the atomic core (a possible way to lift this approximation 
is considered in Section |^), one can present the Green function via a complete set of 
the Volkov wave functions 0p(r,t) that account for the electron dressing by the laser 
field 



(2vr)= 



0p(r,t) = exp\i 
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(5) 



The multiphoton nature of the problem makes the phases of the integrand in to vary 
rapidly with t, t'. This circumstance allows one to use the saddle-point approximation 
for integrations over the time variables. After the accurate integration over the momenta 
p in (0), see details in ||, 0, it is possible to show that (|]), (H), (H) lead to the following 
convenient presentation of the HHG amplitude 
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Each term in is written as a product of two amplitudes of physical, fully accomplished 
and observable processes; no "off-shell" entities appear. Therefore Eqs. (P), (|^) clearly 
demonstrate the stepwise character of the process. The first step is described by the first 
factor Amno(Km) which is an amplitude of physical ATI process when after absorption 
of m laser photons the active electron acquires a translational momentum p = K^. In 
the Keldysh-type approximation this amplitude can be evaluated using the saddle-point 
technique discussed in Section ^. The subscript /iq in A^fMo specifies the contribution 
of one particular saddle point in t' integration, see more details in Section |[ The other 
factor, B]\frn^io(^m), IS a combined amplitude of the second and third steps, which are 
the propagation and laser assisted recombination (PLAR). It can be further factorised 
into the propagation factor 1 / Rmfj.o describing the second step and the amplitude of the 
third step, which is the laser assisted recombination (LAR), CNm(Km)- 
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Rm/10 is merely an approximate expression for the distance passed by the active electron 
in course of its laser- induced wiggling motion 
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The amplitude CNm{p) of the physical LAR process describes the laser assisted recom- 
bination, i.e. transition of an electron with momentum p from the continuum to the 



bound state. Since the continuum state is laser-dressed, the recombining electron can 
emit the A^-th harmonic photon, gaining necessary extra energy from the laser field. 
More details on LAR amplitude are given in Section ^. 

The summation in formula (|^) runs over a number of photons m absorbed on the 
first step of the process when the active electron is released. In other words one can say 
that the energy conservation constraint selects the discrete set of ATI channels in the 
laser-dressed continuum. These channels serve as intermediate states for the three-step 
HG process. In a given channel the electron has translational momentum with the 
absolute value defined by 



Km = V2 {muj -Up + Ea) . (10) 

Here Up = F"^/ (4ci;^) is the well-known ponderomotive potential. ATI plays a role of 
the first stage of HHG process only if the electron momentum has specific direction, 
namely K^, is directed along F. This ensures eventual electron return to the core that 
makes the final step, LAR, possible as discussed in detail in Refs. P, §|, |3- 
The observable HG rates TZn are expressed via the amplitudes as 
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Q = Nuj is the frequency of emitted harmonic, c is the velocity of light. 

Eq. is the major result of this section. It presents the amplitude of the "com- 
plicated" HG process in terms of the amplitudes of "elementary" processes which are 
the ionization and LAR. Conceptual significance of this result is based on the fact that 
it supports the three-step interpretation of the HG process which, in turn, stems from 
the physical picture of the atomic antenna discussed in Section |I|. Eq. (||) is also very 
convenient for numerical applications, since a knowledge on sufficiently simple elemen- 
tary amplitudes enables one to calculate accurately the HG process, as discussed below 
in Section |. 



3 Effective channels 

The summation over m in Eq. (^) has a clear physical interpretation. After the first 
step (ionization) the released electron can be found in any ATI channel before colliding 
with the parent atom. After the third step all intermediate channels result in the same 
final electron state. Therefore the contributions of intermediate ATI level interfere, 
as shows the summation over m in Eq. (|^). The interference has several prominent 
manifestations including the cutoff of the HHG rates for high N and their oscillations 
in the plateau domain, as discussed in Section |^. The number of the ATI channels 
which give significant contribution to the amplitude is usually large 

(5m > 1 . (12) 

This circumstance does not pose a problem for numerical applications based on Eq. (^), 
but may in some situations obscure the qualitative analyses. To overcome this difficulty 
it is desirable to carry out summation over m in (^) in an analytical form. This is the 
major task of this Section which follows the approach of Ref. |]10| . 



Let us verify first that the summation over m in (P) can be replaced by integration 
over the related continuum variable 

dN = Y. '^Nm - / dNm dm • (13) 

To prove the reliability of the approximation based on (0) we use the Poisson sum- 
mation formula which allows one to present the amplitude of HHG in the following 
form 

J2d%m = J2 dm exp(-27rijm) d^^ . (14) 

m j 

The summation index j in (|T^) may be looked at as a variable which is conjugate 
to the m variable. Since muj refers to the spectrum, 27ij/uj should be identified as a 
time variable or, more accurately, as a number of periods of time that elapsed between 
the ejection of the electron from an atom and its return back. The wave function of 
the released electron spreads in space the more the larger j is. Therefore one should 
anticipate that the most important contribution to (|14D is given by the major term 
j = 0. Same conclusion can be drawn from Eq. (|13|). The estimate shows that the 
spectral variable m covers a wide region which causes the time variable j to be localized. 
This discussion demonstrates that one can safely take the leading j = term in (|14D 



thus supporting (O). 



An advantage of integration over the variable m in (|T^) is that it can be carried out 
in closed analytical form using the saddle-point method. To specify this statement let 
us return to Eqs. (||), (0), (^. One can deduce from them that the major dependence 
of the phase of the integrand of @ on integration variables is associated with the factor 
exp{iS) where S is the classical action 



S = S(t,t',m) = f dr 
Jt' 
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Since we know from Eq. (|^) that the momentum variable arises in the final formulas 
as Km, we can use this momentum in ([ISD instead of the integration variable p that 
originates from (^. It should be noted that Eq. (|^) where arises was derived using 
accurate integration over all momenta p; hence the substitution p — > in (plSj) is not 
an additional approximation. 

Eq. (|^) was obtained using the saddle points approximation for integrations over 
the time variables t, t'. Positions of m-dependent saddle points tm, t'm governed by 
the equations 

^5(t,t',m) = 0, (16) 
|^5(t,t',m) = 0, (17) 

in which m is considered as an integer labeling the physical ATI channel. As was 
mentioned above, Eq. ([T3| ) opens a possibility for integration over m that can also be 
carried out using the saddle point approximation. The positions of corresponding saddle 
points are governed by the following equation 

^5(t,t',m) = , (18) 
am 



where one can write the partial derivative over m because the m-dependence of tm and 
does not contribute due to (|16|), (|1^. Eqs. (|1^, (p^, (|1^) define two instants of 
time tm, t'^ at which, respectively, the electron emerges from an atom and returns back 
to it, as well as the number of quanta m = rngfr absorbed in course of the ionization. 
All these three variables are, generally speaking, the complex-valued functions of the 
frequency Vt = Nuj of the generated harmonic. For a given there can be several 
solutions of Eqs. (|16D, (0), (||). 

Integrating over m in Eq. ([l3|) by the saddle method gives the following representa- 



tion for the HHG amplitude 

= 2 E ( 7^ ) dNm.. > (19) 




which is the major result of this Section. Comparing (|19|) with one finds, along with 
clear similarities, several distinctions. The most substantial of them originates from 
different physical meaning of the summation index in these formulas. In Eq. (P) it is 
an integer labeling channels in the physical ATI spectrum. In contrast, in formula (p!9[) 
summation runs over complex- valued m^^ set. It is natural to say that mefj are labels of 
effective channels. In order to find the amplitude dj^^.^^ one can use representation (|^ 
for dj^^ and continue the amplitudes Am^io(i^m) and BNm^ioiJ^m) into the complex-m 
plane which can be done if they are known sufficiently well. One more, though less 



important difference, is presented by an additional square root factor in (19) that arises 
from integration over m and depends on the second derivative S'^ of the action (|T3p 
over m. 

The advantage of effective channel representation stems from the fact that only small 
number of effective channels (actually one or two) contributes, whereas the number of 
essential real ATI channels is quite large, as discussed above. Bearing this in mind it 
is worthwhile to illustrate variation of effective channels labels rrics with A^ by solving 
numerically set of equations ([17|) , (|16D, (|18D for some particular case [[l^]. Fig. 1 shows 



two important solutions mes{N) that move along trajectories in the complex-m plane 
as A^ varies. The overall picture comprises a characteristic cross-like pattern. For small 
A^ the trajectories are close to the real axis and are well separated. They approach each 
other as A^ increases and almost "collide" at some particular critical value A^ = Nc- 
For larger N, N > N^, the trajectories start to move almost perpendicular to the real- 
m axis and rapidly acquire large imaginary parts. It can be demonstrated that large 
Im(meff) lead to suppression of the HHG process. Therefore the critical Nc marks the 
beginning of the cutoff region for HHG. As shown in Ref. |T^, approximate analytical 



solution of Eqs. {^i7\) , (|T6D, (|T8D shows that the critical value Nc is equal to 

NcU =\Ea\+ 3.1731 Up , (20) 
in agreement with the well known result of Refs. P, ^ ^, p2 



In order to find simple physical interpretation for the effective channels, let us note 
that Eqs. ([T6|) , (p^Tf), ([T8|) may be considered as classical equations of motion in the laser 
field. They define the classical trajectories along which the electron first goes away from 
the atom, and then returns back accumulating during this motion energy from the laser 
field that is necessary for HG. This physical picture agrees with the atomic antenna 
concept discussed in Section |1]. It also comes in line with the Corcum model |0] based 
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Figure 1: Trajectories of the effective channel labels mes{N) in the complex m-plane for 
laser frequency to = 0. 0043a. u., intensity I = 10^^ W/cm^ and varying harmonic order N 
for HHG by H~ ion. Positions of two effective channel labels for odd integer are denoted 
respectively by circles and diamonds. 



entirely on the classical trajectories. Even more close relation can be found with the 
approach of Lewenstein et al , which uses the saddle point method to integrate over 
the momenta p in formula (^. 

Eqs. (H) and (|19D provide two different ways to describe the atomic antenna concept, 
either in terms of real physical channels in the intermediate ATI spectrum, or by using 
the effective channels for the intermediate state. Each approach has its advantages 
which can be beneficial for different aspects of HHG problem. Importantly, the above 
discussion ensures identity of the two formulas since (|T9|) was derived directly from (^). 



4 Photoionization and recombination 

Eq. shows that the process of HHG is intimately related to ionization and LAR. 
This fact makes the latter processes very interesting from the perspective of the atomic 
antenna concept, in addition to their well known importance as the basic events in the 
laser-matter interaction. This Section describes the recent progress in the theory of 
these two "elementary" phenomena. Consider first the multiphoton ionization. Adopt- 
ing the Keldysh-type approach which neglects the field of the atomic core in the final 
state one can present the amplitude of the ionization in the following form 

A„(p™) = ^ j\<PpJt) I V{t) I Mt)) , (21) 

where (pp(t) = 0p(r, t) is the Volkov wave function (|) and pm satisfies the energy 
conservation law 

Ea + mLU = ^pl + Up. (22) 



Fast variation of the phase of the integrand in (plD allows one to use the saddle-point 
approach for integration over the time t. This approximation, first proposed by Keldysh 



TT|, was developed in detail in Refs. [12, 13, 14]. Using this scheme one presents the 
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Figure 2: Detachment of ion in bichromatic field with the frequencies to = 0.0043 
a.u. and 3u) and intensities Ii = lO^'^W/cm^ and I2 = lO^W/cm^ respectively. Differential 
detachment rate (in units 10~^a.u.) as a function of the electron emission angle 9 is shown 
for various values of the field phase difference if as indicated in the plots. Open symbols show 
results of calculations of [|l9| (in the = i^vr plot the open circles show the results for ip = ^tt 



and open triangles these for ip = — ^vr). Solid curves show results of the adiabatic theory |18] 



(which coincide for ip = ^ir and (p = — ^vr). Left - the first ATD peak, corresponding to 
absorption of n = 8 photons of frequency uj. Right - the third ATD peak, corresponding to 
absorption of n = 10 photons of frequency iv. 



photoionization amplitude as 

^m(Pm) = ^^m/x(Pm) ; (23) 

where summation runs over essential saddle points labeled by subscript fi. Note that 
in Eq. (^, relating the ionization problem with HG, it is sufficient to take into account 
the amplitude which arises from only one saddle point (of two) labeled as /io; another 
saddle point gives the same contribution that produces a factor 2. 

The technique described above was refined in Refs. [^, |l6l. In the pioneer publi- 



cations |Tl], |T2|, |T3|, |n| the momentum p of the ionized electron was treated as a small 
quantity, and all essential functions were expanded in powers of p/ k. This approxima- 
tion restricted the applicability of the approach, since for high channels in ATI spectrum 
the momentum is not small. Refs. [^, demonstrated that the technique can be mod- 
ified to include large electron momenta as well. Importantly, this modification retains 
simplicity and clear physical nature of the Keldysh approach. 

One can anticipate that the Keldysh-type methods should produce reliable results 
for photodetachment of negative ions, where the detached electron is infiuenced neg- 



ligibly by the Coulomb field of the residual atomic particle. Due to this reason the 
photodetachment was in the focus of attention of Refs. |T^, |T6| which compared the re- 
sults of improved Keldysh approximation with a variety of numerical and experimental 
data available for negative ions. Refs. [|r^, ^ continued this study and extended it to 



the case of the two-color laser field. Detailed description of all results obtained in these 
works would bring us too far away from the main topic of this paper. However, it is 
important to mention that the overall accuracy of the modified Keldysh approximation 
proves be very high. It closely reproduces results of other, much more sophisticated 
methods for total probabilities of detachment as well as for spectral and angular distri- 
butions of photoelectrons both for the weak and strong field regimes (i.e., for any value 
of the Keldysh parameter 7 = ku/F). An example of photodetachment of the H~ in 
bichromatic laser field with the frequencies u = 0.0043 a.u. and Su shown in Fig. 2 
illustrates this point. 

Let us now turn our attention to the other relevant problem, laser-assisted photo 
recombination (LAR). Consider the electron-atom impact in a laser field which results 
in creation of a negative ion and HG. Since the system can acquire energy from the laser 
field due to absorption of several laser quanta, the emitted harmonics should exhibit 
the equidistant spectral distribution. Strange enough, this important process was not 



studied theoretically until two recent almost simultaneous publications [^, pi 
The amplitude of the photorecombination can be written as 

Cm{p) = ^ Hut) I exp{int) U |0p(t)) . (24) 
1 Jo 

Here U = e-r describes the potential which is responsible for the harmonic production, 
and Q is the energy of the generated harmonic which satisfies the energy conservation 
constraint 

n = Y + Up+\Ea\+muj , (25) 

in which m is the number of laser photons absorbed during recombination. The wave 
function of the electron in the continuum in Eq. (p^ can be described by the Volkov 
wave function, similarly to the Keldysh-type approach to photoionization. Comparing 
the amplitudes (0) and (|2T|) one observes their close similarity. This fact allows one 
to develop the theory on the basis of formula (^4|) along the lines described above for 
the ionization problem. In particular, one can use the saddle-point approximation for 
integration over the time variable in (|2^). 

This approach, suggested in Ref. |2^, is supplemented in the cited paper by sev- 
eral numerical examples. One of them, shown in Fig. 3, depicts the cross section of 
recombination on hydrogen atom in a laser field with u = 0.0043 a.u. and the intensity 
/ = 10^^ W/cm^ versus the energy of the emitted high harmonic. Since this is the first 
work in the field we could not compare our results with other calculations. Bearing 
in mind that the recombination process has similarities with the ionization problem, 
where similar approach works well one can expect reliable results in the recombination 
problem as well. 
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Fi gure 3 1 Cross section (Jmip) for laser- assisted recombination of the electron with the 
energy E^i = 10 eV into the bound state in H~ ion (20] . The results are shown for the laser 
field with the frequency uj = 0.0043 a.u. and the intensity / = 10^^ W/cm^. The symbols are 
joined by lines to help the eye. 

5 Quantitative illustrations for HHG 

To illustrate the applicability of the two methods, the factorization technique and the 
method of effective channels, consider an example of HHG by hydrogen negative ion 
in a laser field with oj = 0.0043 a.u. . Using the factorization procedure one needs 
first to calculate the amplitudes of ATI and LAR, which can be done by the technique 
discussed in Sections ^. After that employing (|D one finds the amplitude of HHG, and 



from (|TT|) the HHG rates. The results are presented in Fig. 4. 

It is important to note a strong interference of contributions coming from differ- 
ent intermediate channels m. In the plateau region it is responsible for an oscillatory 
pattern, and becomes even more important in the cutoff region, where a contribution 
of any single channel drastically exceeds the results of an accurate summation over 
a large number of channels. In order to apply the technique based on the effective 
channels one needs to calculate the "elementary" amplitudes for the complex- valued 
number rries of quanta absorbed on the first step of three-step process. This can be 
done by the approach described in Section |^, because it relies on analytical methods 
for calculation of these amplitudes that remain valid for a complex-valued m. Taking 
the effective channels (presented in Fig. 1 for I = 10^^ W/cm^), calculating for them 
the "elementary" amplitudes and applying formula (|l^), we find the rates presented 
in Fig. 4. If one takes into account in the summation ( P^Uj) over mes a single effective 
channel, shown by diamonds in Fig. 1, then the cutoff region is nicely described, as well 
as the overall pattern in the plateau domain. However such one-saddle-point approxi- 
mation does not reproduce structures in the iV-dependence of HG rates. Taking into 
account the two effective channels, shown by diamonds and circles in Fig. 1 improves 
the results for the rates in the plateau domain by producing appropriate structures in 
the A^- dependence. Remarkably, this two-saddle-point calculation gives correct posi- 
tions of minima and maxima in the rate A^-dependence, albeit the magnitudes of the 
rate variation is reproduced somewhat worse; for instance the depth of the minimum 
at A^ = 17 is quite strongly overestimated. 




Figure 4: Harmonic generation rates (|TT|) (in sec^^) for ion in the laser field with the 
frequency uj = 0.0043 a.u. and various values of intensity / as indicated in the plots. Closed 
circles - results obtained by Becker et al |22], open circles - calculations ^ based on factor- 
ization procedure and performing numerical summation (^) over contributions of different ATI 
channels, open diamonds - calculations of Ref. |2C] based on the effective channels approach 
( [T9I) with a single effective channel mefr(A^) taken into account (namely, the effective channels 
shown by diamonds in Fig. 1); open squares - same but taking into account two saddle points 
mc{N) (namely, the effective channels shown by diamonds and circles in Fig. 1). The symbols 
are joined by lines to help the eye. 



Fig. 4 shows good agreement of the approaches based on the factorization technique 
and on the effective channels, which both are in accord with the results of Ref. p2| . 
This agreement holds both for multiphoton regime (left part of Fig. 4), as well as in 
the tunneling regime (right part of Fig. 4). 



6 Above-Threshold Ionization in high channels 

The methods discussed in Sections ||, ^ can be applied to a number of other multiphoton 
problems. To illustrate this point consider the factorization technique for ATI. The 
ionization amplitude An{pn) in the Keldysh-type approximation, considered in (pl|), 
neglects interaction of the released electron with the core. Let A^^pn) be a correction 
that takes this interaction into account. In this notation the total amplitude of ATI 
with absorption of n photons is Aj;°*(p„) = An{pn) + ^n\Pn)- Using the approach of 
Ref. we find the following relation 

= EE E Am{crK^)^f{Pn,crKj , (26) 

which presents the sufficiently complicated correction A^\pn) in terms of two "elemen- 
tary" amplitudes Am(Km) and /^(p, K). The later one describes the electron-atom 
impact in the laser field in the Born approximation. This scattering can be named 
quasielastic, since the atom remains in the same state, but the electron momentum is 
changed K ^ p both in direction and absolute value. Summation over /i in formula 
(p6|) refiects the fact that the electron emission into the continuum takes place at two 




Figure 5: Differential above-tlireshold detachment rates Wn{0)- Left - the rates versus the 
number n of absorbed photons for fixed value of electron ejection angle 9. Right - same as 
functions of for fixed n. An H~ ion is irradiated by the laser wave with frequency lo = 0.0043 
a.u. and intensity / = 5 x 10^'^W/cm^. The minimum number of photons necessary for 
ionization is rimin = 11- Curves - present calculation with account of rescattering; closed circles 
- Keldysh-type approximation of [15|. The arrow indicates the plateau cutoff as predicted by 
classical theory 



moments of time (/i = 1, 2) per laser period. The electron return to the atom is 
ensured only if the electron momentum is parallel (for one value of /i) or antiparallel 
(for another /i) to the field, see details in Ref. This fact is taken into account by a 
summation index cr = ±1 in (|26|). 

Eq. ([26| ) is obtained via application of the factorization technique to ATI. It has 
a transparent physical meaning. Ionization with absorption of n photons needs that 
first m quanta are absorbed by an atom removing the electron from an atom into the 
continuum state with momentum ±Km. The collision of this electron with the atom 
(often referred to as rescattering) results in absorption of additional n — m quanta 
and transition of the electron into a state with momentum p„. All intermediate ATI 
channels labeled by index m contribute coherently. This physical picture of ATI agrees 
with the atomic antenna concept, as was first discussed in [Q. 

We applied (^) to calculation of Above Threshold Detachment (ATD) from H" ion. 
Fig. 5 shows the results which clearly indicate that the contribution of the process (pGf ) 
to the angular distributions of ATD spectra is dominating for higher ATD channels 
while for low channels rescattering effects are small . 



7 Eikonal approach to the Coulomb field 

The Keldysh-type approximation for laser induced ionization, as described in Section 
^, discards interaction of receding photoelectrons with the atomic core. The most 
significant part of this interaction arises due to the Coulomb field of the core. That is 
why a number of numerical applications considered above is carried out for negative ions, 
where the core is neutral. However, for the multiphoton processes with neutral atoms, 
the Coulomb forces between the active electron and the residual positive ion become 
operative. This interaction was taken into account first in Ref. [|1^ in the quasistatic 



limit of the small Keldysh parameter = nuo / F <^ 1. A very simple relation was 
obtained between the photoionization rate wc for the electron bound by a potential 
with the Coulomb field produced by a core charge Z and its counterpart Wsr for the 
electron with the same binding energy but bound by short range forces: 

Wc= [ — \ Wsr ■ (27) 

This result means that for conventional conditions the presence of the Coulomb field 
enhances the rates by several orders of magnitude. It is remarkable that the relation 
(p7|) holds in fact for arbitrary value of the Keldysh parameter, both in the multiphoton 
and tunneling region, as was established in a more elaborate theory by Perelomov and 
Popov IjT^. This result agrees well with experimental data for the total rates ||23|| . 

The theory of Perelomov and Popov is restricted to ejection of low-energy electrons. 
It is usually anticipated that these electrons give principal contribution to the total 
rates summed over all ATI channels as well as over ejection angles. The current exper- 
iments, however, are able to select an individual ATI channel even for a large number 
of absorbed quanta. Both energy and angular distributions of these electrons mani- 
fest some fascinating features which are the object of interest in modern experiment 
and theory. This fact prompts to develop a theory which, reproducing the Perelomov- 
Popov results for low electron energy, could also describe high-energy photoelectrons. 
The antenna-type phenomena considered in this paper provide an additional inspira- 
tion for this study. The ionization amplitude in Eqs. (^, (|26|) should be summed over 
the number m of quanta absorbed that should be large enough. This makes the elec- 
tron momentum in the intermediate ATI channels to be also large. Therefore we need 
to know how the Coulomb field affects the amplitudes for large momenta. In order 
to address this issue we develop the eikonal approach for this problem. It presents a 
simplified version of the semiclassical approximation, that assumes that the Coulomb 
field does not produce significant distortions of classical trajectories describing electron 
propagation in the continuum (i.e., the electron wiggling motion in the laser field). The 
Coulomb field comes into the picture through its contribution to the action / {Z/r)dt 
calculated along the classical trajectories discussed. Within the semiclassical approx- 
imation the action plays a role of the wave function phase. It is important that the 
Coulomb field can produce large contribution to the action while its distortion of the 
trajectories can remain unsubstantial Calculating the action, one is able to construct 
the semiclassical wave function for the released electron, and find with its help the ion- 
ization amplitude. If we totally neglect in this scheme the Coulomb field putting Z = 0, 
the eikonal wave function simplifies to be the Volkov function, and we return to the 
Keldysh-type approximation. An important verification of our eikonal approach pro- 
vides the limit of low photoelectron energy where we reproduce the Perelomov-Popov 



result (p7|). Fig. 6 presents a quantitative example which illustrates importance of the 



Coulomb interaction. 



8 Conclusions 



Atomic antenna gives a clear physical idea how the complicated multiphoton processes 
are operative. Absorption of large number of quanta from the laser field needs that 



1=1 

7 




10 20 30 40 5t> eo 70 &0 90 



0304 



0,003 



11=13. Eik 

11=13. PP 



i!i 0,003 



0,001 



0,000 




10 20 30 40 50 60 70 BO 90 

& (its) 




10 30 30 40 50 80 70 eo 00 
9 (d£g) 



O.OOi 



O.tKJS - 



D.OOE 



0,001 



O.tKW 



■ v. 

: \ 


n=14, 

n=14, 


Eik 
PP 


'. 1 

\ 






\ 



















ID HO aO 40 &D 80 70 AG 90 
9 (dfig) 



Figure 6: Angle-resolved photoionization rates Wn{0) for four lowest open ATI channels 
labeled by a number of absorbed photons n (n = 11 corresponds to the lowest open ATI 
channel). Hydrogen atom is illuminated by the laser wave with frequency oj = 2 eV and 
intensity I = lO^^W/cm^; solid curves - eikonal theory, dashed curves - the Perelomov-Popov 
theory H (|^. 

one of the atomic electrons is at first released from an atom. After that, propagating 
in the core vicinity, the electron accumulates high energy from the laser field and, 
returning to the atomic core, transfers this energy into other channels such as HHG, 
ATI or others. Importantly, this physical picture is implemented in a simple and reliable 
formalism. We discussed two convenient ways to present the theory. One of them, 
called the factorization technique, is presented by Eqs. (^), ( p6D for the cases of HHG 
and ATI. In this approach the amplitude of a complicated process is expressed via the 
physical amplitudes of more simple, "elementary" processes. Another scheme, called 
the effective channels method, is based on Eq. (|19D for HHG. The effective channels are 
closely related to the classical trajectories, that makes them convenient for qualitative, 
as well as numerical studies. We demonstrated an equivalence of the two approaches. 

Applications of both formalisms need calculation of the "elementary" amplitudes. 
This can be achieved by using the modified Keldysh-type approach which very ac- 
curately reproduces data for photodetachment, and, hopefully, for recombination and 
electron-atom scattering in the laser field as well. The Coulomb field of the atomic core 
can be taken into account within eikonal approximation. Reliability of the theoretical 
approaches is demonstrated by quantitative applications. 
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